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O . A Fermi gas of atoms with resonant interactions is predicted to obey univer- 
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sal hydrodynamics, where the shear viscosity and other transport coefficients 



j^ ! are universal functions of the density and temperature. At low temperatures, 

CN . 

v^ . the viscosity has a universal quantum scale hn where n is the density, while at 

(^ ! high temperatures the natural scale is p|//i^ where pr is the thermal momen- 

O. 

O . tum. We employ breathing mode damping to measure the shear viscosity at 

low temperature. At high temperature T, we employ anisotropic expansion of 
^ '. the cloud to find the viscosity, which exhibits precise T^/^ scaling. In both ex- 

periments, universal hydrodynamic equations including friction and heating 
are used to extract the viscosity. We estimate the ratio of the shear viscosity to 
the entropy density and compare to that of a perfect fluid. 

Ultracold strongly interacting Fermi gases are of broad interest, as they provide a tunable 
tabletop paradigm for strongly interacting systems, ranging from high temperature supercon- 
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ductors to nuclear matter. First observed in 2002, quantum degenerate, strongly interacting 
Fermi gases are being widely studied (I7]|2]|2]E1)- To obtain strong interactions (characterized by 
a divergent s-wave scattering length), a bias magnetic field is used to tune the gas to a broad col- 
lisional (Feshbach) resonance, where the range of the collision potential is small compared to 
the interparticle spacing. In this so-called unitary regime, the properties of the gas are universal 
functions of the density n and temperature T. The universal behavior of the equilibrium ther- 
modynamic properties has been studied in detail ([5] |^ HID 13 ESI CU), whereas the measurement 
of universal transport coefficients presents new challenges. 

The measurement of the viscosity is of particular interest in the context of a recent conjec- 
ture, derived using string theory methods, which defines a perfect normal fluid (|72)) . An example 
of a nearly perfect fluid is the quark-gluon plasma produced in gold ion collisions, which ex- 
hibits almost perfect frictionless flow and is thought to be a good approximation to the state of 
matter that existed microseconds after the Big Bang (|73l) . The conjecture states that the ratio 
of the shear viscosity i] to the entropy density s has a universal minimum, r]/s > h/^A-KkB). 
This ratio is experimentally accessible in a trapped unitary Fermi gas, where the entropy has 
been measured both globally ^^ and locally UOWll) and the viscosity can be determined from 
hydrodynamic experiments ([74] [22] E^ 123), so that the predicted minimum ratio can be directly 
compared to that from Fermi gas experiments (|76l[77l) . 

In a Fermi gas, the rf/ s ratio for the normal fluid is expected to reach a minimum just above 
the superfluid transition temperature (I76l) . This can be understood using dimensional analysis. 
Shear viscosity has units of momentum/area. For a unitary gas, the natural momentum is the 
relative momentum hk of a colliding pair of particles, whereas the natural area is the resonant s- 
wave collision cross section (18), An /k"^. Thus, r] oc fik^. At temperatures well below the Fermi 
temperature at which degeneracy occurs, the Fermi momentum sets the scale so /c ~ 1/L, 
where L is the interparticle spacing. Then r] oc fi/L^ and r] oc hn. For a normal fluid above the 



critical temperature, the scale of entropy density s ~ nkB, so r]/s ~ h/kB- For much higher 
temperatures above the Fermi temperature, one expects that hk is comparable to the thermal 
momentum p^ = \/2nikBT, giving the scale r] oc p^/h^ oc T^/'^ /h^ . 

To properly measure the shear viscosity with high precision over a wide temperature range, 
we use universal hydrodynamic equations, which contain both the friction force and the heating 
rate, to extract the viscosity from two experiments, one for each of two temperature ranges. 
For measurement at high temperatures, we observe the expansion dynamics of a unitary Fermi 
gas after release from a deep optical trap and demonstrate the predicted universal T^/^ tem- 
perature scaling. For measurement at low temperatures, we employ the damping rate of the 
radial breathing mode, using the raw cloud profiles from our previous work (|79|) . The smooth 
joining(discontinuity) of the data from the two measurement methods when heating is in- 
cluded(excluded) (|2Ql) . demonstrates the importance of including the heating as well as the 
friction force in the universal hydrodynamic analysis. 

The experiments employ a 50-50 mixture of the two lowest hyperfine states of ^Li, which 
is magnetically tuned to a broad Feshbach resonance and cooled by evaporation in the optical 
trap. The initial energy per particle E is measured from the trapped cloud profile (|2Q|) . 

In the high temperature regime, the total energy of the gas E is larger than 2 Ep, well above 
the critical energy Ec < 0.8 Ep for the superfluid transition (l9l[70l[77]). In this case, the density 
profile is well fit by a Gaussian, n{x, y, z, t) = no(t) exp(— x^/a^ — y'^/o-y — z^ /a1), where 
ai{t) is a time dependent width and no(t) = N/{'n^/'^a,j.ayaz) is the central density and N is the 
total number of atoms. 

The aspect ratio a'x{t)/az{t) is measured as a function of time after release to characterize 
the hydrodynamics, for different energies E between 2.3 Ep and 4.6 Ep, Fig. [T] We also take 
expansion data at one low energy point E = 0.6 Ep, where the viscosity is small compared to 
that obtained at higher temperatures and the density profile is approximately a zero temperature 



Thomas-Fenni distribution. The black curve shows the fit for zero viscosity and no free param- 
eters. To obtain a high signal to background ratio, we measure the aspect ratio only up to 1.4. 
For comparison, the green dashed curve shows the prediction for a ballistic gas. 

We determine the shear viscosity r] by using a hydrodynamic description of the velocity field 
v(x, t) in terms of the scalar pressure and the shear viscosity pressure tensor, 

(9, + vV)t;. = /. + E^^^^^, (1) 
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where f = —VP/n is the force per particle arising from the scalar pressure P and m is the atom 
mass. For a unitary gas, the bulk viscosity is predicted to vanish in the normal fluid (|27]|22]) . so 
we do not include it in the analysis for the expansion. The second term on the right describes the 
friction forces arising from the shear viscosity, where ctjj = dvi/dxj + dvj/dxi — 25ijV ■ v/3 
is symmetric and traceless. 

For a unitary gas, the evolution equation for the pressure takes a simple form, since P = 
2£/3 (|25ll2?l) . where £ is the local energy density (sum of the kinetic and interaction energy). 
Then, energy conservation and Eq.[I]implies (9i + v- V + 5V- v/3)P = 2g/3. Here, the heating 
rate per unit volume q = r] (y1jj2 arises from friction due to the relative motion of neighboring 
volume elements. To express this in terms of the force per particle, /j, we differentiate this 
equation for P with respect to Xj, and use the continuity equation for the density to obtain 

Force balance in the trapping potential t/trap(x), just before release of the cloud, determines the 
initial condition /j(0) = diUtrapi^)- 

These hydrodynamic equations include both the force and the heating arising from viscos- 
ity. The solution is greatly simplified when the cloud is released from a deep, nearly harmonic 
trapping potential Utrap, as /i(0) is then linear in the spatial coordinate. If we neglect viscosity. 



the force per particle and hence the velocity field remain linear functions of the spatial coordi- 
nates as the cloud expands. Thus 9i(V ■ v) = and the pressure P does not appear in Eq. |2l 
By numerical integration (l25l) . we find that non-linearities in the velocity field are very small 
even if the viscosity is not zero, because dissipative forces tend to restore a linear flow profile. 
Hence, the evolution equations [U and [2l are only weakly dependent on the precise initial spatial 
profile of P and independent of the detailed thermodynamic properties. 

We therefore assume that the velocity field is exactly linear in the spatial coordinates. We 
take fi = ai{t)xi and ai{t) = bi{t)ai{0), i.e., the density changes by a scale transformation (|26l) . 
where current conservation then requires Vi = Xi bi(t)/bi{t). 

Li general, the viscosity takes the universal form r] = a{9) Tin, where 6 is the local reduced 
temperature and 77 — )• in the low density region of the cloud (|2 71 1201) . Using the measured trap 
frequencies, and eqs. [T]and|2l the aspect ratio data are fit to determine the trap averaged vis- 
cosity parameter, a = [l/Nfi) J (i^x?7(x, t), which arises naturally, independent of the spatial 
profile of 77. Since 9 has a zero convective derivative everywhere (in the zeroth order adiabatic 
approximation) and the number of atoms in a volume element is conserved along a stream tube, 
a is a constant that can be compared to predictions for the trapped cloud before release. 

As shown in Fig. [H the expansion data are very well fit over the range of energies studied, 
using a as the only free parameter. We find that the friction force produces a curvature that 
matches the aspect ratio versus time data, while the indirect effect of heating is significant in 
increasing the outward force, which increases the fitted a by a factor of ~ 2, compared to that 
obtained when heating is omitted (|20|) . 

For measurements at low temperatures, where the viscosity is small, we determine a from 
the damping rate of the radial breathing mode (fTPJ) . For the breathing mode, the cloud radii 
change by a scale transformation of the form 6j = 1 + e^, with Cj << 1, and the heating rate in 
eq. [2]is oc ef, which is negligible. Hence, the force per particle evolves adiabatically. Adding 



the trapping force to eq.[Tl one obtains the damping rate 1/r = ha/{3m{x'^)) (|20l [ 

The fitted viscosity coefficients a for the entire energy range are shown in Fig. 2, which 
can be used to test predictions (|29]li3l22])- Despite the large values of a at the higher energies, 
the viscosity causes only a moderate perturbation to the adiabatic expansion, as shown by the 
expansion data and the fits in Fig. [H The breathing mode data and expansion data smoothly 
join, provided that the heating rate is included in the analysis. In contrast, omitting the heating 
rate produces a discontinuity between the high and low temperature viscosity data (|2Q|) . The 
agreement between these very different measurements when heating is included shows that 
hydrodynamics in the universal regime is well described by eqs.[T]and|2l 

To test the prediction of the T^/^ temperature scaling in the high temperature regime, we 
assume that r] relaxes to the equilibrium value in the center of the trap, but vanishes in the low 
density region so that a is well defined. This behavior is predicted by kinetic theory (|27|) . We 
expect that a ~ ao where r]o = ao^^o is the viscosity at the trap center before release. At high 
temperatures (|75l) . 

ao = as/2 9o, (3) 

where 03/2 is a universal coefficient. As 6 has a zero convective derivative everywhere (in the 
zeroth order adiabatic approximation), 6^0 at the trap center has a zero time derivative and ao is 
therefore constant as is a. 

The inset in Fig. [2] shows the high temperature (expansion) data for a versus the initial 
reduced temperature at the trap center, 6*0. Here, 6q = TQ/Tp{nQ) = {TQ/Tpi){nj/no)'^^^. The 
local Fermi temperature Tpino) = h^ (Sn'^no)'^/^ / {2mkB) and Tpi = Ep/kB = Tpini) is the 
ideal gas Fermi temperature at the trap center, nj is the ideal gas central density for a zero 
temperature Thomas-Fermi distribution. We use {ni/nof'/'^ = 4(o-^/cr|.^)/7r^/'^ and obtain the 
initial Tq/Tfi from the cloud profile (|20|) . 

The excellent fit of Eq.[3]to the data, inset Fig.|2l demonstrates that at high temperature, the 



viscosity coefficient very well obeys the 9q scaling, in agreement with predictions (|75l) . We 
note that Eq.[3]predicts that a^ scales nearly as E^, because 6*0 oc Tq/hq oc E"^ . This explains 
the factor of ~ 10 increase in the viscosity coefficients as the initial energy is increased from 
E = 2.?,EFioE = A.QEF. 

A precise comparison between the viscosity data and theory requires calculation of the trap- 
average a from the local shear viscosity, where the relation is tightly constrained by the ob- 
served T^/^ scaling. Our simple approximation a ~ ao, yields «3/2 = 3.4(0.03), where 0.03 is 
the statistical error from the fit. A better estimate based on a relaxation model (|32| ) shows that 
a = 1.3 ao at high T, yielding 03/2 = 2.6. At sufficiently high temperature, the mean free path 
becomes longer than the interparticle spacing, since the unitary collision cross section decreases 
with increasing energy. In this limit, a two-body Boltzmann equation description of the viscos- 
ity is valid. For a Fermi gas in a 50-50 mixture of two spin states, a variational calculation (|75l) 
yields 03/2 = 45 7r^/^/(64-\/2) = 2.77, in reasonable agreement with the fitted values. 

Finally, Fig. [3] shows an estimate of the ratio of rj/s = ahn/s = {h/kB)a/{s/nkB) ^ 
{h/kB)a/S, where S is the average entropy per particle of the trapped gas in units of ks- We 
obtain S in the low temperature regime from Ref. ^, which joins smoothly to the second virial 
coefficient approximation for S in the high temperature regime ([20l) . The inset shows the low 
temperature behavior, which is about five times the string theory limit (red dashed line) near the 
critical energy ([9|) Ec/Ep = 0.7 — 0.8 (|20|) . We note also that the apparent decrease of the r]/s 
ratio as the energy approaches the ground state (J9l) 0.48 Ep does not require that the local ratio 
— )■ as T — 7- 0, since contributions from the cloud edges significantly increase 5* compared to 
the local s at the center. 
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Figure 1: Anisotropic expansion. (A) Cloud absorption images for 0.2, .3, 0.6, 0.9, 1.2 ms 
expansion time, E = 2.3 Ep; (B) Aspect ratio versus time. The expansion rate decreases at 
higher energy as the viscosity increases. Solid curves: Hydrodynamic theory with the viscosity 
as the fit parameter. Error bars denote statistical fluctuations in the aspect ratio. 
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Figure 2: Trap averaged viscosity coefficient a = f cPx r]/{hN) versus initial energy per atom. 
Blue circles: Breathing mode measurements; Red squares: Anisotropic expansion measure- 
ments. Bars denote statistical error arising from the uncertainty in E and the cloud dimensions. 
Inset: a versus reduced temperature 6^0 at the trap center prior to release of the cloud. The 
blue curve shows the fit oq = 0:3/2 Oq ' demonstrating the predicted universal high temperature 
scaling. Bars denote statistical error arising from the uncertainty in 9q and a. A 3% systematic 
uncertainty in Ep and 7% in ^0 arises from the systematic uncertainty in the absolute atom 
number (l20l) . 



12 



10-1 



6- 



05 



4- 



0- 



0.8 
0.6 
0.4 
0.2 
0.0 



ii 



^ 



"ffr 



_ ■ 



0.6 0.7 0.8 0.9 1.0 1.1 



m 



m )t 



f* 



^* 



«»w 



i' 



t« 



*l 



0gg 



I* 



* 



e 



s 



T r 

2 3 

E/E. 



-i 1 

4 5 



Figure 3: Estimated ratio of the shear viscosity to the entropy density. Blue circles: Breathing 
mode measurements; Red squares: Anisotropic expansion measurements; Inset: Red dashed 
line denotes the string theory limit. Bars denote statistical error arising from the uncertainty in 

E, a, and S ^. 
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